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ABSTRACT 

We present 10 new gamma-ray burst (GRB) redshifts and another five redshift limits based on host 
galaxy spectroscopy obtained as part of a large program conducted at the Very Large Telescope (VLT). 
The redshifts span the range 0.345 < z < 2.54. Three of our measurements revise incorrect values 
from the literature. The homogeneous host sample researched here consists of 69 hosts that originally 
had a redshift completeness of 55% (with 38 out of 69 hosts having redshifts considered secure). 
Our project, including VLT/X-shooter observations reported elsewhere, increases this fraction to 77% 
(53/69), making the survey the most comprehensive in terms of redshift completeness of any sample 
to the full Swift depth, analyzed to date. We present the cumulative redshift distribution and derive 
a conservative, yet small, associated uncertainty. We constrain the fraction of Swift GRBs at high 
redshift to a maximum of 14% (5%) for z > 6 (z > 7). The mean redshift of the host sample is 
assessed to be (z) > 2.2, with the 10 new redshifts reducing it significantly. Using this more complete 
sample, we confirm previous findings that the GRB rate at high redshift (z > 3) appears to be in 
excess of predictions based on assumptions that it should follow conventional determinations of the 
star formation history of the universe, combined with an estimate of its likely metallicity dependence. 
This suggests that either star formation at high redshifts has been significantly underestimated, for 
example due to a dominant contribution from faint, undetected galaxies, or that GRB production is 
enhanced in the conditions of early star formation, beyond that usually ascribed to lower metallicity. 
Subject headings: dust, extinction — galaxies: distances and redshifts — gamma rays: bursts 



1. INTRODUCTION 

Determining the statistical properties of gamma-ray 
bursts (GRBs) has long been compromised by inhomoge- 
neous selecti on and a bias again st optically dark bursts. 
With Swift (Gehr els et al.ll2004l ) it has become possible 
to construct much more uniform samples, and to target 
the host galaxies even of optically faint bursts via X-Ray 
Telescope (XRT) localizations, for which redshifts could 
not be determined from the afterglows. 

We have been securing GRB host galaxy information 
for a homogeneous sample of 69 Swift GRBs with a 
large program at the Very Large Telescope (VLT). The 
first observations of The Optically Unbiased GRB Host 
(TOUGH) sample were obtained on 2006 February 23 
with the survey concluding on 2008 August 29. The im- 
mediate goals are to determine the host luminosity func- 
tion (LF), study the effects of reddening, determine the 
fraction of Lya emitters in the hosts, and obtain red- 
shifts for targets without a reported one. The sample 
has been carefully selected and obeys strict and well- 
defined criteria. To optimize the survey, we focused on 
systems with the best observability, which also have the 
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best available infor mation. The survey design and cata- 
logs are presented in lHjorth et al.l (l2012f). the fundamen- 
tal properties of t he hosts in [M alesani et all (12012). the 
Lya emission in IMilvang- Jensen et al. (2012), and new 
VLT/X-shooter redshifts in[Kruhler et al.l (|20 121 ). 

This paper presents the first TOUGH campaign for 
missing redshifts. We attempted spectroscopic observa- 
tions of most TOUGH host candidates with R < 25 mag 
that did not have a reported reliable redshift. A to- 
tal of 19 candidates were spectroscopically observed for 
an aggregate of nearly 30 hr with the aim of acquir- 
ing redshift information. In addition, we include a 
new redshift measurement for GRB 06090 8 reported i n 
IMilvang- Jensen et al.l (|2012l see also lFvnbo et al.ll2QQ9h . 
The details of our observations and reductions are de- 
scribed in the next section. New redshifts and redshift 
limits for each observed host are presented in ^3l We 
model the redshift distribution and compare it to the ob- 
served one in §4\ Finally, the implications of our observa- 
tional efforts are discussed in $5) We adopt a cosmology 
where the Hubble paramet er is Hp = 70.4 km s -1 Mpc -1 , 
tt m = 0.27 and tt A = 0.73 dJarosik et all 1201 11 1. 

The reduced data from this work will be available from 
ESO0 and from the TOUGH Web sited 

2. OBSERVATIONS AND DATA REDUCTION 

The host spectroscopic observations were carried 
out between 2006 May 30 and 2008 August 29. 
The FORS2 instrument was used as well as FORS1 
(jAppenzeller et all [1998) before and after the blue 

6 http://archive.eso.org 

7 http : //www . dark-cosmology . dk/TOUGH 
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CCD upgrade. Four different grism setups were 
mainly used: 300V, 300V+GG375, 600z+OG590, and 
6OORI+GG4350 We requested a seeing constraint of l'/2 
for all observations apart from GRB 060923C for which 
0'/8 was required ( §3.13p . For the majority of the ob- 
servations, this goal was accomplished. More detailed 
information is listed in Table [TJ 

The data reduction was performed following standard 
techniques for bias and flat-field corrections. The indi- 
vidual spe ctra were then cosm ic ray cleaned using the 
method of Ivan Dokkuml (|2001l ). If the host trace was 
clearly visible in individual images, the spectra were op- 
timally extracted for each two-dimensional image sep- 
arately. Otherwise, the two-dimensional images were 
aligned, averaged, and the host spectrum extracted if 
detected. 

The wavelength calibration was applied using a He- 
NeAr lamp spectrum obtained in the morning after the 
science observations. The root mean square scatter 
around the wavelength calibration fit was roug hly 0.2 A 
for the 300V grism and less than 0.1 A for the higher 
resolution grisms. The individual wavelength-calibrated 
spectra were averaged, and the corresponding Poisson 
error spectra (calculated by the IRAF/apall task) were 
quadratically averaged. Finally, flux calibration was ap- 
plied by using standard star observations carried out 
every night a host was observed. We note, however, 
that the transparency constraint was set to "thin cir- 
rus" implying that absolute flux calibration should be 
interpreted with caution. 

3. NEW REDSHIFTS AND LIMITS 

A total of 20 host systems were observed (including 
GRB 060908), of which only 10 had a reported optical 
afterglow (OA) or near-infrared afterglow (NIRA). In 
only two cases (GRBs 060805 A and 070808) there is more 
than one galaxy detected within or on the border of the 
XRT error circle. For the other eight systems without 
an OA/NIRA the host identity is nearly unambiguous 
with the probability of cha nce pro jection being f airly low 
(frequently less than 5%; Ma lesani et al.l 12012). Below 
we discuss each system in detail and justify our redshift 
measurements and limits. When no spectral features are 
detected we set redshift limits in the following way: 

— If the continuum is visible we estimate, using the 
error spectrum, the lowest wavelength at which it is sig- 
nificantly detected (at a significance level of around 3a). 
The lack of Lya forest lines is then used to set a firm 
redshift upper limit. 

— To get a redshift lower limit we note that GRB hosts 
are in general star-forming galaxies and thus display the 
usual emission lines, such as [O II] A3727, H/J, [O III] 
A5 007, and Ra; s ee e.g. the GHostS Web sit e F j Table 3 
in iSavaglio et al.l (|2009h . and iLevesque et al.l (|2010f ). In 
particular, the non-detection of [O II] A3727 in our 300V 
spectra immediately sets a redshift limit of z > 0.9. This 
value is based on 7200 A, since above it there are strong 
skylines which leave significant residuals in the reduced 
spectra. 

8 See http://www.eso.org/sci/facilities/paranal/ 
instrument s/fors/doc/ for more details. 

9 http://www.grbhosts.org/ 



— If neither emission lines nor a continuum are de- 
tected we follow lGrupe et a l. (2007) and assign a redshift 
upper limit to bursts with excess (above Galactic) X-ray 
absorbing column densityQ calculated at z = 0, above 
an equivalent hydrogen colum n density of 2 x 10 2 * cm -2 
(including the uncertainty)!^ iGrupe et al.l (j2007l ) used 
an upper limit of z = 2, but we will be slightly more con- 
servative and assign an up per limit of z = 3 .5 to these 
bursts (see Equation (1) in lGrupe et al.ll2QQ7h . The cor- 
responding re st-frame column densi ty is 10 23 cm -2 . In 
the sample of iCampana et al.l (j2010l . see their Figure 2), 
no burst exceeds this number, lending support to our 
choice of z < 3.5. This approach was utilized for six 
events in the whole TOUGH sample. Among those, five 
have an i?-band host detection, implying that their red- 
shifts have to be z < 6. 

3.1. GRB 0507 UB (No OA/NIRA) 

There is a single faint source detected inside the XRT 
error circle, with a bright source close to the southern 
edge of the error circle. Our spectroscopic observations 
show the brighter one to be an M star. Neither emission 
lines nor a continuum are detected from the fainter ob- 
ject and hence no redshift information can currently be 
obtained for this host candidate. However, we conclude 
that z < 3.5 based on the excess column density detected 
in the X-ray spectrum. 

3.2. GRB 050822 (No OA/NIRA) 

In the 300V spectrum, the continuum is clearly de- 
tected in the 3900-7200 A region corresponding to the 
redshift range 0.9 < z < 2.2 (left panel of Fig. QJ. There 
is an emission line candidate (7a) at 9071 A (right panel 
of Fig. CQ) that is most likely the [O II] A3727 doublet at 
z = 1.434. In support of this, we note that the spec- 
tral extent of the line is significantly broader than the 
extent of emission lines visible in other traces in the two- 
dimensional spectrum. 

At this redshift, our spectra would not cover H/3, [O III] 
A5007, and Ha. Other potential line identifications (of 
9071 A) such as H/3, [O III] A5007, and Ha are rejected 
on account that we would expect to see additional and 
roughly similarly bright lines in the wavelength range 
covered. 

3.3. GRB 0509 15 A (NIRA) 

The host galaxy of the dark (Jakobsson et al. 2004) 
GRB 050915 A was previ o usly r eported in lOvaldsen et al.l 
(|2007T ) and lPerlev et al.l (|2QQ9f ). In 300V the continuum 
is clearly detected in the 4300-7200 A region without 
unambiguous features, corresponding to 0.9 < z < 2.5 
(Fig.©. 

In the upper panel of Fig. [2j the continuum of another 
brighter galaxy on the slit is visible (z = 0.444)13 It 

10 We use the afterglow-only late spectrum (photon counting 
mode) obtained from http://www.swift.ac.uk/xrt_spectra/ 

11 For compar i son, o ne of the sample selection criteria, as denned 
in Hiorth et al. (2012), is that the Galactic extinction is Ay < 
0.5 mag, which corresponds to a hydrogen column dens ity of 0.9 x 
10 21 cm" 2 according to the relation derived by Predehl & Schmitt 
(p95h . 

12 This value was erroneously reported as the GRB 050915A host 
redshift in Jako bsson et al . ( 2009, 2 0lTallbllcI) . 
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Figure 1. One- and two-dimensional spectra (left: 300V; right: 600z) of the GRB 050822 host. There are no emission lines identified in 
300V. Telluric features and skyline residuals are marked with 0, whereas the error spectrum is plotted as a dotted line. The 600z emission 
line candidate is marked with a circle in the right panel. Hashed regions mark areas strongly affected by skylines. 
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Figure 2. One- and two-dimensional spectra (300V) of the 
GRB050915A host. In the upper panel, it is the bottom trace 
between the two horizontal dashed lines. The Lya absorption fea- 
ture is marked with a vertical line, whereas telluric features and 
skyline residuals are marked with 0. The error spectrum is plotted 
as a dotted line. 
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Figure 3. One- and two-dimensional spectra (300V) of the 
GRB 051001 host. In the upper panel, it is the bottom trace be- 
tween the two horizontal dashed lines. The Lyoi absorption feature 
is marked with a vertical line, whereas telluric features and skyline 
residuals are marked with 0. The error spectrum is plotted as a 
dotted line. 

clearly contains flux down to around 3600 A. Therefore, a 
spectral break must be present in the GRB 050915 A host 
continuum, which we interpret as the Lya break. More 
specifically, the ID spectrum shows a flux drop around 
4300 A corresponding to z « 2.54. Recent observations 
by VLT/X-shooter confirm this interpretation where we 
detect [O II] A3727, H/3, a nd [O HI] AA4959, 50 07 in emis- 
sion at a similar redshift (Kriihler et al. 2012f). 

3.4. GRB 051001 (No OA/NIRA) 

The continuum is extremely faint in the 600z spectrum 
and there are no signs of any emission lines. In 300V 
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Figure 4. One- and two-dimensional spectra (600RI) of the 
GRB 051006 host. Emission lines are marked with vertical lines, 
whereas telluric features and skyline residuals are marked with 0. 
The error spectrum is plotted as a dotted line. 

the continuum is clearly detected in the 4200-7200 A 
region without unambiguous features, corresponding to 
0.9 < z < 2.5 (Fig. p. Thus, GRB 051001 can safely 
be ruled out as a hig h-redshift burst as suggested by 
iSalvaterra et al.l ()2QQ7T ). 

The ID spectrum shows a flux drop around 4170 A 
which we interpret as the Lya break at z ~ 2.43. Re- 
cent observations by VLT/X-shooter confirm this inter- 
pretation where we detect H/3 , [O HI] A5007, and Ha in 
emission at a similar redshift (Kriihle r et al.ll2012[ ). 

3.5. GRB 051006 (No OA/NIRA) 

The 600RI spectrum clearly shows a strong and broad 
emission line which we identify as the [O II] A3727 dou- 
blet at z = 1.059. At that redshift, we also detect a much 
weaker line, [Ne III] A3869 (Fig. gj). The stronger line is 
unlikely to be H/3, [O III] A5007, or Ha since additional 
and roughly similarly bright lines would be expected in 
the wavelength range covered. 

3.6. GRB051117B (No OA/NIRA) 

Our 300V spectrum clearly shows a strong emission 
line in addition to a few absorption features (Fig. [5]). 
These features are consistent with being [O II] A3727, 
H10, H9, H8, He and Ca II K and H at a common redshift 
of z = 0.481. We note that Re and Ca II H are blended 
with a strong skyline residual. 

3.7. GRB 060306 (No OA/NIRA) 

Neither emission lines nor a continuum were detected 
in the 600z spectrum and hence no redshift information 
could be obtained. However, there is a faint detection of 
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Figure 5. One- and two-dimensional spectra (300V) of the 
GRB051117B host. Line features are marked with vertical lines, 
whereas telluric features and skyline residuals are marked with 0. 
The error spectrum is plotted as a dotted line. 
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Figure 6. One- and two-dimensional spectra (300V) of the 
GRB 060306 host. Telluric features and skyline residuals are 
marked with 0, whereas the error spectrum is plotted as a dot- 
ted line. The diamonds show the rebinned (100 A) spectrum. 
The vertical dashed lin e indicates the location of Lya if z = 3.5 
(|Salvaterra et al.ll20T3 ). 

the continuum in the 300V spectrum in the approximate 
wavelength range 4300-6800 A (Fig. E| indicating 0.8 < 
z<2.5. 

In particular, as apparent from the rebinned spectrum 
(diamonds in Fig. 0), we do not detect any flux break 
corresponding to Lya at z = 3.5. This redshift was based 
on the detection of a single e mission line interpreted as 
[O II] (jSalvaterra et al.ll2012l ). We note that this lack 
of a break is different from what we detect from z > 2 
hosts presented in this paper (GRBs 050915 A, 051001, 
and 070129). We tentatively suggest that the redshift 
must be either z = 2.35 or z = 1.55 if the emission line 
is interpreted as [O III] or Ha, respectively. 

3.8. GRB 060719 (NIRA) 

The continuum is very faint in 600z but slightly 
brighter in 300V (Fig. [7j). Based on the continuum in 
the 3600-7200 A region we infer 0.9 < z < 2.0. Indeed, 
recent VLT/X-shooter observations have detected emis- 
sion lines (strong Ha and an indication of [O II] A3727) 
outside of our grism wav elength range, cons istent with 
the FORS redshift limit (jKruhler et alJl2()Tl . 

3.9. GRB 06080 5 A (No OA /NIRA) 

There are two host galaxy candidates within the XRT 
error circle, a bright object (A) at the southwestern edge 
and a s econd, fainter source (B) slightly northe ast of the 
center (jPerlev et aLll2QQ9b iMalesani et al.ll2012l ). Our slit 
covered both objects. No line features are observed over 
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Figure 7. One- and two-dimensional spectra (300V) of the 
GRB 060719 host. In the upper panel, it is the fainter upper trace. 
Telluric features and skyline residuals are marked with 0, whereas 
the error spectrum is plotted as a dotted line. 
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Figure 8. One- and two-dimensional spectra (300V) of the 
brighter galaxy (object A) within the XRT error circle of 
GRB 060805A. The fainter galaxy (object B) is barely seen in the 
two-dimensional spectrum just above object A. Line features are 
marked with vertical lines, whereas telluric features and skyline 
residuals are marked with 0. The error spectrum is plotted as a 
dotted line. 

the spectral range for object B. There is, however, a faint 
continuum detectable by block averaging the spectrum 
along the dispersion axis. It is detected down to approx- 
imately 4200 A which corresponds to z < 2.5. 

The continuum of object A is detected down to at least 
3600 A whi ch implies z 5 i 2.0 (F ig. [8]), consistent with the 
findings of iPerley et al.l ([2009). The spectrum also dis- 
plays [O II] A3727 and [O III] A5007 at a common redshift 
of z = 0.603. The probability of chance projection for ob- 
jects A and B is 1.6% and 6 .6%, respectively, c alculated 
following the prescription in Bloom et all (|2002f ). Hence, 
the host identification is ambiguous, and we cannot claim 
to have secured the redshift of this burst. 

3.10. GRB 060814 (NIRA) 

The host system of GRB 060814 se ems to consist of 
three different blobs (Malesani 2006) which we mark 
as A, B, and C ( Fig. ©. The location of the NIRA 
([Levan et al.l 120061) is consistent with blob A; we used 
the data from lLevan et al.l ([20061 ) along with our images 
to precisely positio n the afterglow on the host complex. 
iThone et all (|2007f ) reported a redshift of z = 0.84 al- 
though it is not clear which part of the complex was 
covered by their slit. 

Our spectra clearly show emission lines emanating 
from the combined region of blobs B and C (Fig. [TO]) . 
In the 600z spectrum (seeing: l'/3), we detect H/3 and 
[O III] AA4959,5007 at a common redshift of z = 0.841. 
There is an indication of [O II] A3727 in the 300V spec- 
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Figure 9. i?-band image (seeing: 0'/65) of the GRB 060814 field. 
The location of the NIR afterglow is marked with a cross and is 
consistent with blob A. The orientation and extent of the FORS 
slit is also shown; it mostly cover blobs A and B, and partly blob 
C. The X-shooter slit only covers blobs A and B. 
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Figure 10. Two-dimensional spectrum (600z) of the GRB 060814 
host complex. The dashed upper (lower) line indicates the centroid 
of blob A (B+C) which is shown in Fig. [9] Hashed regions mark 
areas strongly affected by skylines. 

trum but residuals from strong sky lines make the iden- 
tification ambiguous with our resolution. 

However, there are no indications of any emission lines 
at the location of blob A. A clear detection of the con- 
tinuum in the 300V spectrum in the approximate range 
3600-7200 A suggests that 0.9 < z < 2.0. In any case, 
the conclusion is that a host redshift of z = 0.84 is 
not supported by our data. Indeed, recent VLT/X- 
shooter observations led to the detection of [O II] A3727, 
[O HI] AA4959.5007 , and Ra in emission at z = 1.92 
(|Kriihler et al.l 120121 ). clearly visible in blobs A and B 
(the slit did not cover blob C as shown in Fig. [9|)F^I At 
this redshift, the most prominent emission lines are out- 
side of our FORS grisms wavelength range. The nec- 
essary consequence of these observations taken together 
is that the host complex consists of blobs A and B at 
z = 1.92, while blob C is a chance alignment at z = 0.84. 

3.11. GRB 060908 (OA) 

A Lya emission line is detected from the host galaxy 
at a re dshift of z = 1.884 (part of t h e TOUGH Lyq cam - 
paign: iMilvang- Jensen et aLll2Q12h - iFvnbo et all (|2009l ) 
utilized these findings and confirmed the redshift by de- 
tecting the C IV AA1548, 1550 doublet in the OA spec- 
trum, ruling out the tentative redshift of z = 2.43 re- 

13 Our analysis i s consistent with what is reported by 
I Sal vat err a et al.l (|2012T ) who used the same dataset. 
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Figure 11. Two-dimensional spectrum (300V) of the 
GRB 060923 A host. It has been block averaged along the 
dispersion axis to make the host continuum visible. It is the 
bottom trace between the two horizontal dashed lines. 

ported bv lRol et all (|2006l ). 

3.12. GRB 060923 A (NIR A) 

The ext remely red af t erglow of GRB 060923A was dis- 
cussed in iTanvir et al.l ([2008). Their later-time optical 
imaging and spectroscopy revealed a faint galaxy coinci- 
dent with the GRB position which implied a moderate 
redshift and therefore that dust is the likely cause of the 
very red afterglow color. 

Our 300V spectrum reveals a very faint continuum 
with a possible weak emission line candidate at around 
6330 A (3cr detection). In order to verify this we obtained 
a higher resolution 1200R spectrum. There is again a 
weak detection of the continuum but no sign of an emis- 
sion line at the aforementioned wavelength. In Fig. [11] 
the continuum is detected dow n to approximately 4600 A 
which corresponds to z < 2.8 (jTanvir et al.l 2008). 

3.13. GRB060923C (NIRA) 

The GRB 060923C NIRA was located only around l'/5 
away from an R = 20.9 mag bright star. Detecting the 
host (R — 25.5 mag) spectrum, even under good seeing 
conditions, is therefore quite a challenge. As can be seen 
in the upper panel of Fig. [12j the stellar spectrum over- 
whelms any signs of the host continuum. However, there 
is an indication of an emission line at the expected loca- 
tion of the host trace at approximately 6942 A. In order 
to get a better view of the situation we subtracted the 
star continuum using a two-dimensional deconvolution 
technique. 

The deconvolution method applied to th e spectrum is 
presented in detail in iCourbin et al.l ((2000), resulting in 
a spatially deconvolved spectrum with the stellar spec- 
trum separated from the spectrum of the host galaxy. 
This spatial separation of the pointlike component from 
the diffuse one relies on an accurate building of the point- 
spread function and on the hypothesis that the host con- 
tains no significant structure narrower than the fixed and 
finit e resolution of the deconvolved spectrum. Follow- 
ing iLetawe et al.l (|2QQ8[ ) , we have chosen to subtract the 
stellar spectrum, as derived from the deconvolution pro- 
cess, from the original spectrum. This avoids smoothing 
of the diffuse component mandatory for a proper sep- 
aration between the star and host. Thus, we obtain a 
two-dimensional host galaxy spectrum with a spatial res- 
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Figure 12. Original (upper panel) and deconvolved (lower panel) 
two-dimensional spectrum (300V) of the GRB 060923C host. The 
emission line candidate is indicated with a circle. 

olution constant along the slit and equal to that of the 
original data 

The results are shown in the lower panel of Fig. [12] 
The emission line candidate is still visible. If real, the 
identification could only be with Lya or [O II] A3727 
in which case the redshift would be z = 4.71 or z = 
0.86, respectively. We have considered other potential 
line identifications such as H/3, [O III] A5007, or Ha all 
of which are rejected on account that we would expect 
to see additional and roughly similarly brig ht lines in the 
wavele ngth range covered. We note that the lGrupe et al.l 
(2007) method implies that z < 3.5 so that z = 4.71 is 
excluded. 

3.14. GRB 061021 (OA) 

Our 300V host spectrum shows [O II] A3727 and [O III] 
AA4959, 5007 at a common reds hift of z = 0.345 (Fi g. [13]). 
Using this redshift information, iFvnbo et al.l ([2009) iden- 
tified the Mg II AA2796, 2803 doublet in the very blue end 
of the afterglow spectrum. 

3.15. GRB 070103 (No OA/NIRA) 

The host continuum detection down to 7500 A in our 
600z spectrum only allows us to set a modest limit of 
z < 5. 2 (Fig. [T4l). This is less constraining than us- 
ing the iGrupe et al.l (120071) method which implies that 
z < 3.5 (|Fvnbo et al.l 120091 ). There are no convincing 
emission lines present. However, recent VLT/X-shooter 
observations have detected [O III] AA4959, 5007 in emis- 
sion outside of our grism wavelength ra nge, consistent 
with the aforementioned redshift limits (Krii hler et al.l 

EoIH). 

3.16. GRB 070129 (OA) 




4000 5000 6000 o 7000 8000 

Observed wavelength [A] 



Figure 13. One- and two-dimensional spectra (300V) of the 
GRB 061021 host. Emission lines are marked with vertical lines, 
whereas telluric features and skyline residuals are marked with 0. 
The error spectrum is plotted as a dotted line. 




7500 8000 8500 9000 9500 



Observed wavelength [A] 

Figure 14. One- and two-dimensional spectra (600z) of the 
GRB 070103 host. Telluric features and skyline residuals are 
marked with 0, whereas the error spectrum is plotted as a dot- 
ted line. 




4000 5000 6000 o 7000 8000 

Observed wavelength [A] 



Figure 15. One- and two-dimensional spectra (300V) of the 
GRB 070129 host (lower trace). The Lya absorption feature is 
marked with a vertical line, whereas telluric features and skyline 
residuals are marked with 0. The error spectrum is plotted as a 
dotted line. 

Neither emission lines nor a continuum were detected 
in the 600z spectrum and hence no redshift information 
could be obtained. However, there is a clear detection of 
the continuum in the 300V spectrum between approxi- 
mately 4100 and 7200 A, suggesting that 0.9 < z < 2.4. 
Thus, we can immediately rule GRB 070129 out as a 
high-r edshift candidate as suggested by iSalvaterra et al.l 
(|2007h . 

In the upper panel of Fig.[15j the continuum of another 
similarly bright object on the slit is visible. It clearly 
contains flux down to around 3900 A. Therefore, a spec- 
tral break is present in the GRB 070129 host continuum, 
which we interpret as the Lya break. More specifically, 
the ID spectrum shows a flux drop around 4070 A cor- 
responding to z « 2.35. Recent observations by VLT/X- 
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Figure 16. One- and two-dimensional spectra (300V) of the 
GRB070419B host. Telluric features and skyline residuals are 
marked with 0, whereas the error spectrum is plotted as a dot- 
ted line. 

shooter confirm this interpretation where we detect [O II] 
A3727, [O IH l AA4959, 5007, and R a in emission at a sim- 
ilar redshift (jKruhler et aLll2012[ ). 

3.17. GRB 070306 (NIRA) 

The highly extinguished afterglow of GRB 070306 
and the prope r ties o f its host galaxy were discussed in 
Uaunsen et al.l ([2008). The OA low-resolution spectrum 
displayed a single emission line which was tentatively in- 
terpreted as [O II] A3727. In an attempt to verify this 
interpretation by resolving the doublet we observed the 
host galaxy with the higher-resolution 1 028z grism as 
part of TOUGH. As shown in Figure 4 in Jaunsen et al. 
(2008) the doublet was successfully resolved and a red- 
shift of z = 1.496 confirmed. 

3.18. GRB 070328 (No OA/NIRA) 

A very faint host continuum detection down to 7500 A 
in our 600z spectrum only allows us to set a modest 
limit of z < 5 .2. T his is less constraining than using the 
iGrupe et al.l ([20071) method which implies that z < 3.5 
(jFvnbo et al.l 120091 ) . There are no convincing emission 
lines present. 

3.19. GRB070419B (OA) 

The continuum is detected in the 300V spectrum but 
without any significant emission lines visible (Fig. [T6j) . 
There is an unambiguous detection of the continuum 
between approximately 3900 and 7200 A indicating that 
0-9 ^ z < 2.2. Indeed, recent VLT/X-shooter observa- 
tions have detected [O III] AA4959, 5007 and Ha in emis- 
sion outside of our grism wa velength range, cons istent 
with the FORS redshift limit (jKruhler et al.ll2012h . 

3.20. GRB 070808 (No OA/NIRA) 

There is a faint (R = 26.7) object consistent with the 
X-ray position, as well as two brighter ones just outside 
the X-ray error circle. We obtai ned a 600RI spectrum 
of the brightest source (object A: Malesa ni et al.ll2012f ). 
It clearly shows two emission lines which we identify as 
the [O II] A3727 doublet and H/3 at a common redshift of 
z = 0.681 (Fig.[I7J). The probability of chance projection 
for this object is around 4%-5%, c alculated following the 
prescription in lBloom et al.l ([2002) , moreover it definitely 
lies outside the UVOT-enhanced XRT error circle. The 
fainter object is consistent with the XRT error circle, 
but it also has a large chance superposition probability 




5000 5500 6000 6500 7000 o 7500 8000 
Observed wavelength [A] 



Figure 17. One- and two-dimensional spectra (600RI) of object 
A located close to the NE border of the GRB 070808 XRT error 
circle. Emission lines are marked with vertical lines, whereas tel- 
luric features and skyline residuals are marked with 0. The error 
spectrum is plotted as a dotted line. 

(25%). Hence, we cannot claim to have secured the red- 
shift of this burst. However, we conclude that z < 3.5 
based on the excess column density detected in the X-ray 
spectrum. 

3.21. Summary 

We report the redshift s of 10 GRB host galaxies whose 
average redshift is (z) = 1.59, significantly lower than 
the overall Swift GRB mean redshift. The low value is 
most likely the result of targeting the brightest galax- 
ies in the sample (R < 25 mag). GRB 061021 is among 
the closest "classical" Swift long-duration burst detected 
so far. Furthermore, we have estimated redshift limits 
for an additional five hosts and inferred that three burst 
redshifts reported in the literature are erroneous (GRBs 
060306, 060814, and 060908). The results are listed in 
Table [1] along with the observational details. We have 
also estimated the emission line flux and the equivalent 
width (EW) of most line features detected and detailed 
them in Table [2 

4. THE GRB REDSHIFT DISTRIBUTION 
4.1. Redshift Measurements and Constraints 

Figure [18] shows the cumulative redshift distribution 
of the 53 TOUGH bursts with a measured redshift, in- 
cluding four photometric reds hifts and the n e w VL T/X- 
shooter reds hifts presented in iKriihler et al.l (|2012f ) (see 
Table 3 in IHiorth et~aT1l2012[ ). Also plotted is a conser- 
vative systematic error band (hatched region) containing 
information for all the 69 TOUGH bursts. The shaded 
region represents the likely statistical (la standard error 
of the sample) uncertainty of the measured redshift dis- 
tribution under the assumption that it is a true random 
sample of the overall population. 

The hatched region incorporates the information pre- 
sented in Table [TJ including the redshift limits. In ad- 
dition, the upper boundary is produced by simply plac- 
ing GRBs with only upper limits to their redshifts at 
z = 0. The lower boundary is generated in the follow- 
ing way. GRBs with a detected OA and/or host are 
placed at the maximum redshift they can have, given 
their bluest photometric detection reported in the liter- 
ature or in IHiorth et al.l ([20121 ) . Following the general 
decli ne of the Lya f orest density with decreasing redshift 
(e.g. Songaila 2004), we set two different limits depend- 
ing on the wavelength region probed by the available 
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Table 1 

New Redshifts and Redshift Limits Determined in the Host Galaxy Sample 



GRB 


OA or NIRA? 


Redshift 


^host 


Instrument 


Grism+Filter 


Slit Width 


Seeing 


Spectral Res. 


Exp. Time 










(mag) 






\ ) 


{ ) 


\ A ) 


\ s ) 


050714B 


No 


z 


< 3.5 


25.5 


oFORSl 
nFORSl 


300V 
300V 


1.3 
1.3 


0.9 
0.7 


11.9 
7.3 


4 x 1220 
4 x 1220 


050822 


No 




1.434 


24.2 


oFORSl 
FORS2 


300V+GG375 
600z+OG590 


1.3 
1.0 


0.8 
1.9 


10.6 
3.2 


2 x 1300 
2 x 1275 


050915A 


NIRA 


Z r- 


a 2.54 


24.6 


oFORSl 


300V 


1.3 


1.1 


14.5 


4 x 1295 


051001 


No 


Z r 


a 2.43 


24.4 


nFORSl 
FORS2 


300V 

600z+OG590 


1.3 
1.0 


0.8 
0.7 


8.4 
2.3 


4 x 1360 
2 x 1350 


051006 


No 




1.059 


23.0 


FORS2 


600RI+GG435 


1.0 


1.3 


3.3 


2 x 430 


051117B 


No 




0.481 


21.1 


oFORSl 


300V 


1.3 


1.3 


17.2 


1 x 1000 


060306 


No 


0.8 < z 


< 2.5 


24.1 


oFORSl 
FORS2 


300V+GG375 
600z+OG590 


1.3 
1.0 


1.3 
1.1 


17.2 
3.2 


2 x 1300 
2 x 1275 


060719 


NIRA 


0.9 < z 


< 2.0 


24.6 


nFORSl 
FORS2 


300V 

600z+OG590 


1.3 
1.0 


0.9 
1.0 


9.4 
3.2 


4 x 1360 
4 x 1350 


060805A 


No 


z 


< 2.5 


23.5 (A) 
25.1 (B) 


nFORSl 
FORS2 


300V 

600z+OG590 


1.3 
1.0 


1.3 
1.4 


13.6 
3.2 


4 x 1295 
1 x 1500 


060814 


NIRA 




1.92 


22.9 


nFORSl 
FORS2 


300V 

600z+OG590 


1.3 
1.0 


0.8 
1.3 


8.4 
3.2 


2 x 1295 
2 x 1350 


060908 


OA 




1.884 


25.5 


nFORSl 


600B 


1.3 


0.8 


6.5 


4 x 1345 


060923A 


NIRA 


z 


< 2.8 


26.1 


nFORSl 
FORS2 


300V 

1200R+GG435 


1.3 
1.3 


1.0 
1.0 


10.5 
1.5 


2 x 1280 
2 x 1350 


060923C 


NIRA 


z 


< 3.5 


25.5 


nFORSl 


300V 


1.0 


0.8 


8.4 


4 x 1360 


061021 


OA 




0.345 


24.4 


nFORSl 


300V 


1.3 


0.7 


7.3 


2 x 1280 


070103 


No 


z 


< 3.5 


24.2 


FORS2 


600z+OG590 


1.0 


1.0 


3.2 


4 x 900 


070129 


OA 


Z r 


a 2.35 


24.4 


nFORSl 
FORS2 


300V 

600z+OG590 


1.3 
1.0 


1.4 
1.4 


13.6 
3.2 


4 x 1360 
2 x 1275 


070306 


NIRA 




1.496 


22.9 


FORS2 


1028z+OG590 


1.0 


0.7 


2.5 


2 x 1300 


070328 


No 


z 


< 3.5 


24.4 


FORS2 


600z+OG590 


1.0 


1.0 


3.2 


5 x 1330 


070419B 


OA 


0.9 < z 


< 2.2 


25.2 


nFORSl 


300V 


1.3 


0.9 


9.4 


6 x 1360 


070808 


No 


z 


< 3.5 >26.7 


FORS2 


600RI+GG435 


1.3 


2.2 


4.3 


1 x 715 



No te. — -Rhost i s the -R-band total magnitude (or 3a upper limit) of the host galaxy (before correcting for Galactic extinction) 
from Malesani et al. ( 2012). The GRB 060805A host redshift limit is conservative and refers to the fainter object (B) within the XRT 
error circle. The GRB 060908 host redshift originates fr om the Lyo: pa r t of T OUGH (Milvang- Jensen et al. 2012). The host redshifts 
of GRBs 060719, 070103, and 070419B are reported in IKriihler et al.l (pOll ) and stem from our VLT/X-shooter observations. The 
redshift uncertainty is of the order of ±0.001 for those given to three decimal places. Data obtained with FORS1 before 2007 April 
6 are marked oFORSl. Data obtained after the FORS1 blue CCD upgrade are labeled nFORSl. The wavelength coverage for each 
setup is approximately 3500-8140 A (oFORSl/300V/l{'3), 3500-9640 A (nFORSl/300V/iyO), 3500-8880 A (nFORSl/300V/iy3), 7470- 
10 700 A (FORS2/600z/l"0), 5300-8630 A (FORS2/600RI/1^0), 4950-8250 A (FORS2/600RI/1^3), 7700-9500 A (FORS2/1028z/l"0), 
and 5720-7200 A (FORS2/1200R/1"3). 



observations. At low redshift, the opacity of the Lya 
forest is small, hence, for detections in blue and optical 
filters (V band and blueward), the maximum redshift is 
derived by comparing the reddest wavelength in the filter 
response curve to the Lyman limit. At higher redshift, 
only negligible flux survives blueward of the Lya wave- 
length, t herefore, for redder fil ters, the comparison is to 
Lya (e.g. iDairAglio et alJl2QQ8[ ). If a more constraining 
limit is not available, we can set their maximum redshift 
to z = 3.5 if their X-ray spectra fulfill the criterion de- 
scribed in ^3] For only a single burst, GRB 061004, are 
there no redshift constraints available. It is arbitrarily 
placed at a maximum redshift of z = 10 in Fig. [THJ 

4.2. Modelling 

Various authors have attempted to predict or 
model the GRB redshift distributi on and hence 
comp are it to observations (e.g. iGorosabel et al.l 
20041 : iNataraian et al.l l2005M jakobsson et all I2006F 



Le fe Dermer 



Guetta & Piran 



20071; ISaTyaterra 
MM 



Kistle r et al.l 1200^ „ _ _ 

201(1 Robertson fc Ellis 20ll) 



Tanvir 



Chincarini 
Jakobsson 



2007; 
2007 



Dong fc Lul 120091 : iCampisi et al.l 

but of necessity have 



worked with samples with considerably greater incom- 
pleteness and optical bias than TOUGH now provides. 



During the review process of this work, Salvat erra et al.l 
(|2012[ ) presented the redshift distribution of a complete 
sample of Swift GRBs, which was inspired by criteria 
similar to those applied to TOUGH. Their sample 
is slightly smaller and, by construction, limited to 
gamma-ray bright events. The TOUGH sample extends 
to fainter luminosities (a factor of around six) and as 
such is more suited to test the faint (and potentially 
high-redshift) end of the GRB distribution. Indeed, 
their mean (median) redshift is 1.84 (1.64), significantly 
lower than for the whole TOUGH sample. 

Given the association of lon g GRBs with the 

deaths of massi ve stars (e.g. iHjorth et al.l [2003; 

IWooslev fc Bloom] |2006[ ), it is commonly assumed that 
the GRB rate density fol lows the star-formation rate 
(SFR ) density history (e.g. IPa czvnski 1998 : IWiiers et al.l 
ll998MPorciani fc MadauT 120011 : Uakobsson et al.l 12005: 
iFvnbo et al.l 120081 ), with a possible low- meta llicity 
enhancement (e.g . lLanger fc Normanl 120061 : O 120081 : 
iButler et al.ll201(J ) and/or evolution with redshift in in- 
trinsi c rate or LF (e.g . [Cowar 

3 [2005LE aigne et al. 

2006: IGuetta fc Piranl 120071 : iKistler et al l 120081 : 

iSalvaterra et al.l 120091 ). For the purposes of this 

observational paper, we will only present illustrative 
model fits. 
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Figure 18. Thick solid curve: the cumulative fraction of GRBs as a function of redshift for the 53 Swift bursts in the TOUGH sample 
with a measured redshift ((z) = 2.23). Hatched region: this is a conservative error region showing the systematic error on the thick solid 
curve. Shaded region: statistical region showing the la sampling error band around the thick solid curve. Dotted curve: the expected 
redshift distribution for Swift observable long GRBs using the SFR1 history parameterization (see the main text). Dashed curve: the 
same redshift distribution for the SFR2 history parameterization (see the main text). 



Currently, we have little guidance from theory as to 
plausible functional forms for the GRB LF, and various 
ones have been considered in the above studies. Two 
of the most commo nly adopted forms are either a bro- 
ken power law (e.g. IGuetta et al.l [20051: IGuetta fc Piranl 
I2007t iButler et all 120101 : IWanderman fc Piranl 1201 Of ), or 
a single power la w with an exponentia l cutoff at low 
luminosities (e.g . iNataraian et al.l 120051 : iCampisi et al.l 
l2Q10t iCao et all |2011[ ). Both forms generally under- 
predict the rate of low-luminosity GRBs (as exem- 
plified by GRB 980425), suggesting that such bursts 
may form a separat e population dCobb et al.l 1 2006: 
Soderberg et al.l 120061: LLi ang et al. 2007; Co ward et al.l 
2008: IZitounTet al.l 120081 : iFolev et al.i r2008: Virgili elalj 
2009). In our homogeneous and unbiased TOUGH sam- 



ple, only one event (GRB 060218) is generally regarded as 
a member of this faint population, though we note this 
represents a rate of nearby bursts consistent with that 
predicted from the BATSE sample by Chapman et al. 
(200?]). For the sake of clarity and simplicity, we restrict 
ourselves here to the exponentially cutoff single power- 
law form: 



a (L/L p y 



-L p /L 



(1) 



We assume that the GRB rate follows the SFR 
history, and consider two different SFR history pa- 



rameterizations which we label as follows. SFR1 is 
an update (ID 120081) of the SFR history models of 
Hopkins fc Beacom (2006) to include recent data from 



Bouwens et al. (2008) and iReddy et all (j2QQ8h . com- 



bined with a lo w-metallicity modifi c ation following the 
prescription of Langer & No rmanl ((2006). SFR2 is 
model A from Schmidt (2009) which represents a SFR 
history which remains constant beyond z ~ 3. It 
may, for example, be considered a more extreme low- 
metallicity co rrection to the cosmic SFR, or repres ent 
a correction dKistler et al.l 120091 : iVirgili et all 1201 lft to 
the high-redshift SFR as estimated from flux-limited sur- 
veys (by the integration of galaxy LFs thus obtained) 
due to a large amount of hidden star formation in 
faint, low-mass, and high specific SFR galaxies of the 
type that G RBs tend to be as s ociated with at lower 
redshift (e.g. iFvnbo et all 12001 lLe~Floc'h et al.l |2QQ3 
iChristensen et aLll2QQ4 : iFruchter et al.ll2006ft . 

Modeling is perform ed in the standard manner (e.g. 
IGuetta fc Piranl [2007) to produce log N- log L number 
count distributions for various parameters of the LF, 
which are then fit by y 2 minimization to the observed 
log N-\og L distribution of all Swift bursts with peak pho- 
ton flux > lcm _2 s _1 . We emphasize that the redshift 
distribution is not part of this fitting procedure, but is 
always purely a result. In Fig. [TBI we pl°t the redshift 
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Table 2 

Line Identifications, Observed EWs, and Fluxes 



GRB Feature 


EW obs 

(A) 


Flux (1CT 17 
ergs -1 cm -2 ) 


050822 


[O II 


A3727 


< -23 


2.2 


051006 


O IIJ A3727 
Ne III] A3869 


-95.7 ±4.0 

OQ A _|_ O O 

— Zo.4 zb z.y 


7.5 
2.0 


051117B 


O II] A3727 

H10 

H9 

H8 

Ca II K 

Ga II H and He 


-12.5 ± 1.5 
O.U zb 1. 1 
6.8 ± 1.1 

5.2 ±1.1 


6.6 

— 


060805A 


O II] A3727 
O III] A5007 


-20.0 ±4.0 
-67.0 ±9.0 


2.6 
7.6 


061021 


O II] A3727 
O III] A4959 
O III] A5007 


-49.1 ±9.9 
-24.0 ±9.0 
-51.0 ±9.5 


1.5 
0.8 
1.3 


070306 


[O II 
O II 


A3726 
A3729 




17.0 


070808 


[O II 

H/3 


A3727 


-65.3 ±5.1 
-17.5 ±3.3 


10.9 
4.5 



Note. — Flux measurement errors are of the order of 
20%. When the EW is not given, skylines prohibited a re- 
liable estimate. The fluxes reported for GRB 060805A are 
those from object A. The fluxes reported for GRB 070808 are 
those from object A. The flux in the emission line candidate 
in GRB 060923C is not reported due to the uncertainties in- 
troduced by the deconvolution. 



distributions from our best fitting models in comparison 
to the TOUGH redshift data. The best fit model in- 
corporating the SFR1 parameterization is a good fit to 
the number count distribution with a reduced x 2 = 1-32 
(12 dof) and has parameters L p — 10 49,88 ergs -1 and 
v = 1.82. The resultant redshift distribution is a good 
fit to the TOUGH redshift distribution up to z ~ 1.5, 
but underpredicts bursts at higher redshift. Overall, a 
Kolmogorov-Smirnov (KS) test yields a <0.5% likelihood 
of the observed and model populations being the same. 
The SFR2 parameterization is a slightly better fit to 
the number count distribution with a reduced x 2 = 1-30 



(12 dof) and LF parameters L p = 10 50-23 ergs -1 and 
v = 1.83. The resultant redshift distribution in this case 
is an improved match to the overall observed TOUGH 
distribution with a KS likelihood of 42%. 

At face value, these results seem to imply that GRBs 
follow a cosmic SFR history that is significantly en- 
hanced at high redshift compared to estimates from flux- 
limited surveys. As previously discussed, given what is 
known about GRB hosts, it is entirely feasible that GRBs 
trace star formation at high redshift that would be un- 
detectable by other means. It is of course also possible 
that the simple low-metallicity enhanced SFR parame- 
terization used in the SFR1 model is inadequate, or that 
the LF could have a more complex form and/or evolve 
with redshift. 

An alternative approach to the modeling is to directly 
fit the observed joint peak-luminosity and redshift dis- 
tribution of the sub-sample of TOUGH GRBs with red- 
shifts. Again we take both SFR1 and SFR2, but in 
this case a broken power-law LF. We take a thresh- 
old flux for detection by the Burst Alert Telescope of 
2 x 10 -8 erg cm -2 s -1 , and apply a if-correction based 
on the observed spectral parameters for each burst, to 
bring the luminosity to a fixed restframe bandpass (30- 



300 ke V) . The results are very similar to the conclusions 
obtained above; hence we do not plot the inferred red- 
shift distribution in Fig. [THJ However, we do note that 
for SFR2 the maximum-likelihood solution gives the fol- 
lowing LF parameters: 

<P{L) OC | L -2.00. L > l 52.5 ergs -l W 

5. DISCUSSION AND IMPLICATIONS 

Figure [TBI shows the cumulative redshift distribution of 
the largest homogeneous and unbiased sample of GRBs 
to date. It contains redshift information on all the 69 
TOUGH bursts, including limits. The conservative sys- 
tematic error region (hatched) of the TOUGH sample 
redshift distribution is significantl y smaller than for pre - 
vious samples (e.g. Figure 2 in Jako bsson et al.l 12009). 
This allows the rejection of various model predictions 
(e.g. compare the dashed and dotted curves in Fig. [T8|) . 
Thus, we have been able to confirm previous findings (e.g. 
Kistl er et aLll2QQ9b iVirgili et all 120111 : lR,obertson fc Ellis! 



2012D that the GRB rate at high redshift (z > 3) ap- 
pears to be in excess of predictions based on the as- 
sumption that it should follow conventional determina- 
tions of the star-formation history of the universe, com- 
bined with an estimate of its likely metallicity depen- 
dence (jLanger fc Nor man 2006|). 

It is possible that star formation at high redshifts 
has been significantly underestimated. Even at z ~ 
2 it appears that the galaxy LF ha s a substantially 
steep er faint-end slope than locally (e.g. lReddv fc Steidell 
2009), while recent LF studies in the Hubble Ultra- 
Deep Field have concluded that at z > 7 so-far un- 
detected galaxies are likely to complet ely dominate 
the total star formation activity (e.g. iBouwens et al.l 
120111 : iTanvir et all 1201 2( 1. This picture is supported 
by p revious obse r vations of damped Lyq a b sorbers 



_e 1 g 1 lFvnbo et al.l 119991 : lHaehnelt etHI l2000b IS chave 
2001) and Lya emittin g galaxies (e.g. iFvnbo et stl.ll2003t 



Jakobsson et al.l 120051). as well as by recent simula- 
tion studies ([Choi fc Nagaminell2011f ). Alternatively, it 
could be that GRB production is substantially enhanced 
in the conditions of early star formation, beyond the 
metallicity-dependent rate correction already applied. In 
the long run, large complete samples of GRB redshifts 
should shed light on whether the GRB rate is propor- 
tional to SFR or whether other effects play an important 
role. 

The sampling error and the conservative systematic er- 
ror region are shown separately to clearly illustrate that 
incompleteness dominates the sample, and more is gained 
by reducing the systematics rather than increasing the 
sample size. Using both error regions we can set a con- 
servative limit on the maximum number of Swift bursts 
at z > 6 (z > 7): 14% (5%). This is fully consistent 
with the models which predict between 2 (SFR1) and 
13 (SFR2) bursts per year (all sky) at 6 < z < 9 to 
Swift /BAT limits. 

The average (median) redshift of the 53 TOUGH 
bursts is (z) = 2.23 [z = 2.14), signific antly lower than 
the e arly Swift results indicated (e.g. Jakob sson et al.l 
2006). This difference may simply reflect the compar- 
atively small samples analyzed in that paper, but could 
also be due to an increased success in measuring redshifts 



The Swift GRB Redshift Distribution 



11 



v 




1 2 3 4 5 6 7 



Figure 19. Thick solid curve and shaded region are defined as 
in Fig. 1181 caption. The thin solid curve is the cumulative fraction 
of GRBs as a function of redshift for all long Swift bursts to date 
(2012 April 1: 198 redshifts) with (z) = 2.16. 



available instrume ntation, such as the VLT/X-shooter 
(jVernet et al.ll20lT| ). we can continue to follow this track 
into the future. 

We thank the referee for providing excellent comments 
in a timely manner and Andrew J. Levan for provid- 
ing us with the NIR afterglow image for GRB 060814. 
We also thank Thomas Kriihler, Dan Perley, and Ruben 
Salvaterra for very informative discussions. Based on 
observations collected at the European Southern Obser- 
vatory, Chile, as part of the large program 177.A-0591. 
P.J. acknowledges support by a Project Grant from the 
Icelandic Research Fund. D.M. acknowledges financial 
support from the Instrument Center for Danish Astro- 
physics. The Dark Cosmology Centre is funded by the 
Danish National Research Foundation. This work has 
made use of the University of Hertfordshire Science and 
Technology Research Institute high-performance com- 
puting facility. 



z < 2 using weaker absorption lines in afterglow spectra, 
and via host galaxies. The mean redshift could be as 
low as (z) ~ 1.7 (upper boundary of the hatched region) 
although it is unlikely that the majority of bursts with 
unknown redshifts would be located at very small dis- 
tances. In fact, it is more probable that (z) > 2.20 since 
we have only targeted the brightest galaxies in the sam- 
ple (R < 25 mag) for spectroscopic follow-up. In Fig. [T9j 
we compare the TOUGH redshift distribution to the one 
obtained from all long-duration Swift bursts. The agree- 
ment is good in general but the TOUGH sample in- 
cludes relatively more bursts in the range 1.8 < z < 2.5. 
This can be explained by the fact that Lya is located in 
the blue end of the optical at this redshift range which 
makes it easy to miss a damped Lya absorber. On the 
other hand, emission lines are readily detectable with X- 
shooter at these red shifts. 

As discussed in iHjorth et al.l (j2012[ h one of the 
TOUGH selection criteria (requiring the XRT error circle 
radius to be less than 2") might slightly bias our sample 
against faint, hence potentially high-redshift, bursts. In 
total, only two events were rejected due to a large error 
radius (GRBs 050412 and 061102). In both cases, the 
X-ray afterglows were unusually faint, so this is an esti- 
mate of the fraction^ (2/71) of Swift GRBs that would 
fail the TOUGH positional accuracy selection criterion 
which might well be at high redshift. Further to this, a 
comparable bias is introduced by the exclusion of bursts 
without a detected X-ray afterglow at all (5/180)0 All 
five GRBs are low-significance detections and are thus 
likely to be spurious. In total, a maximum of 5.5% of 
Swift long bursts could be at high redshift and not enter 
the TOUGH sample. 

We have now reached a point in GRB research where a 
single burst rarely elucidates and illuminates our general 
understanding of the field. It is important to focus on 
well-defined samples and population studies, where sys- 
tematics and biases can be minimized. Swift has made 
it possible to build such a sample and thanks to new 

14 The total number of bursts here being the TOUGH 69 along 
with those two with a large error radius. 

15 Swift detected 180 long bursts during the TOUGH time period 
in which the XRT was quickly (< 1 hr) repointed. 
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